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ABSTRACT 

The recently discovered planetary system HD45364 which consists of a Jupiter and Saturn mass planet is very likely 
in a 3:2 mean motion resonance. The standard scenario to form planetary commensurabilities is convergent migration 
of two planets embedded in a protoplanetary disc. When the planets are initially separated by a period ratio larger 
than two, convergent migration will most likely lead to a very stable 2:1 resonance for moderate migration rates. To 
avoid this fate, formation of the planets close enough to prevent this resonance may be proposed. However, such a 
simultaneous formation of the planets within a small annulus, seems to be very unlikely. 

Rapid type III migration of the outer planet crossing the 2:1 resonance is one possible way around this problem. In this 
paper, we investigate this idea in detail. We present an estimate for the required convergent migration rate and confirm 
this with N-body and hydrodynamical simulations. If the dynamical history of the planetary system had a phase of 
rapid inward migration that forms a resonant configuration, we predict that the orbital parameters of the two planets 
are always very similar and hence should show evidence of that. 

We use the orbital parameters from our simulation to calculate a radial velocity curve and compare it to observations. 
Our model can explain the observational data as good as the previously reported fit. The eccentricities of both planets 
are considerably smaller and the libration pattern is different. Within a few years, it will be possible to observe the 
planet-planet interaction directly and thus distinguish between these different dynamical states. 

Key words, planet formation - migration - resonance - N-body simulations - hydrodynamical simulations 



1. Introduction 

Over 400 extrasolar planets have aheady been discov- 
ered (jSchneideil l2009f l and their diversity keeps challeng- 
ing planet formation theory. For example, the recently dis- 
covered multi-planetary system HD45364 raises interesting 
questions about its formation history. 

The planets have masses of mi = 0.1906A/jup and 
1712 = 0.6891Mjup and are orbiting the star at a distance 
of ai = 0.6813 AU and 02 = 0.8972 AU, respectively 
(jCorreia et al.ll2009D . The period ratio is close to 1.5 and 
a stability analysis implies that the planets are deep in- 
side a 3:2 mean motion resonance. The planets have most 
likely formed further out in cooler regions of the proto- 
stellar disc, as water ice which is an important ingredient 
for dust aggregation can only exist beyond the snow line 
which is generally assum ed to be at radii larger than 2 AU 
(|Sasselov fc Lecaill2000t l. 

It is then usually assumed that migration due to planet 
disc interactions has moved the planets closer to the star. 
Although the details of this process are still hotly de- 
bated, the existence of many resonant multiplanetary sys- 
tems supports this idea. During migration the planets can 
get locked into a commensurabilities after which the plan- 
ets migrate together with a constant period ratio. In such 
a re sonance, one or m ore resonant angles are librating (see 
e.g. iLee fc Pealell200lD . 



For the planetary system HD45364 this standard pic- 
ture poses a new problem. Assuming that the planets have 
formed far apart from each other, the outcome for the ob- 
served masses after migration is almost always a 2:1 mean 
motion resonance, not 3:2 as observed. The 2:1 resonance 
that forms is found to be extremely stable. One possible 
way around this is a very rapid convergent migration phase 
that passes quickly through the 2:1 resonance. 

In this paper we explore this idea quantitatively. The 
plan of the paper is as follows. In section [2] we show, us- 
ing N-body simulations, that the system always ends up in 
the 2:1 resonance assuming moderate migration rates. We 
present scenarios which result in formation of a 3:2 reso- 
nance after a rapid migration phase. In section[3]we perform 
hydrodynamic simulations with a variety of disk models in 
order to explore the dependence on the physical setup. In 
section |4] we briefly discuss other formation scenarios. We 
go on to compare the orbital parameters of our simulations 
with the observed radial velocity data in section El 

We find that the orbital parameters observed in our 
simul ations differ from th ose estimated from a statistical 
fit bv lCorreia et al.l (|2009f ). However, our models reproduce 
the observational data very well and we argue that our 
fit has the same lev el of significance as that obtained by 
ICorreia et all (|2009D . 

Future observations will be able to resolve this issue. 
This is the first prediction of orbital parameters for a spe- 
cific extrasolar planetary system derived from planet mi- 
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Fig. 1. Period ratio P2/P1 as a function of time (y-axis) and 
migration timescale of the outer planet Ta,2 (y-axis). The 
migration timescale of the inner planet is Tq^i — 2000 yrs. 
The inner planet is initially placed at ri — lAU. The ec- 
centricity damping is given through K = Ta/T,. = 10. 



gration theory only. The parameter space of orbital config- 
urations produced by planet disc interactions (low eccen- 
tricities, relatively small libration amplitudes) is very small. 
As in case of the GJ876 system, this can provide strong ev- 
idence on how the system formed. Finally, we summarise 
our results in section [6l 



2. Formation of HD45364 

2.1. Convergent migration and resonance capture 

In th e core accretion model (for a review see e.g. iLissauej 
fl99l a solid core is firstly formed by dust aggregation. 
This process is much more efficient if water exists in solid 
form. In the proto-stellar nebula this happens beyond the 
ice line where the temperature is below 150 K at distances 
larger th an a few AU. S ubsequently, after a critical mass is 
attained (|Mizunolll98Cll), the core accre tes a gaseous enve- 
lope from the nebula (|Bodenheimer fc Pollack 1986). Both 
planets in the IID45364 system are interior to the ice line, 
implying that they should have migrated inwards. 

The migration rate depends on many parameters of the 
disc such as surface density and viscosity as well as the mass 
of the planets. The planets are therefore in general expected 
to have different migration rates which leads to the possi- 
bility of convergent migration. In this process the planets 
approach orbital commensurabilities. I f they do this slo wly 
enough, resonance capture may occur (|Goldreichlll9"65l) af- 
ter which they migrate together maintaining a constant pe- 
riod ratio thereafter. 

Studies made by a number of authors have shown that 
when two planets, of either equal mass or with the outer 
one being the more massive, undergo differential convergent 
migration, capture into a mean motion commensurability is 
expected to occur p rovided that the conve rgent migration 
rate is not too fast (ISnellgro ve et al.ll200lh . The observed 
inner and outer planet masses are such that, if (as is com- 
monly assumed for multiplanetary systems of this kind) the 
planets are initially widely enough separated so that their 
period ratio exceeds 2, at low migratio n rates a 2:1 commen- 
surabili ty is expected to form (e.g. iNelson fc Papaloizoul 
1200 2: Kiev et al.ll2004 ) 

fPierens fc Nelsoiil (|2008t ) studied a similar scenario 
where the goal was to resemble the 3:2 resonance between 



Jupiter and Saturn in the early solar system. They also 
found that the 2:1 resonance forms at early stages. However, 
in their case the inner planet had the larger mass whereas 
the planetary system that we are considering has the heav- 
ier planet outside. In this situation the 2:1 resonance can 
be unstable, enabling the formation of a 3:2 resonance 
later on and the migration rate may stall or even reverse 
(Masset fc SneUgrove..2001i) . 

2.2. The 2:1 mean motion resonance 

We found that if two planets with masses of the observed 
system are in a 2:1 mean motion resonance, which has been 
form via convergent migration, this resonance is very sta- 
ble. An important constraint arises, because as indicated 
above, provided the planets start migrating outside any low 
order commensurability, at the slowest migration rates a 2:1 
resonance is expected to form rather than the 3:2 commen- 
surability that is actually observed. 

We can estimate the critical relative migration timescale 
Ta.crit above which a 2:1 commensurability forms from the 
condition that the planets spend at least one libration pe- 
riod while migrating through the resonance. The resonance 
semi- major axis width Aa associated with the 2:1 resonance 
can be estimated from the condition that two thirds of the 
mean motion difference across Aa be equal in magnitude 
to 27r over the libration period. This gives 



Aa = 



"2 



(1) 



where 02 and n2 are the semi major axis and the mean mo- 
tion of the outer planet, respectively. The libration period 
27: /wif c an b e expressed in terms of the orbital parameters 
(see e.g. iGol dreich 1965; Rein & Papaloizou 2009) but is, 
for convenience, here measured numerically. If we assume 
the semi-major axes of the two planets evolve on constant 
(but different) timescales |ai/ai| = Tq^i and |a2/d2| = Ta.2, 
the condition that the resonance width is not crossed within 
a libration period gives 
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to pass through the 2:1 MMR. 

If the planets of the IID45364 system are placed in a 
2:1 resonance with the inner planet being located at 1 AU, 
the libration period 27r/a;;/ is found to be approximately 
75 yrs. Thus, a relative migration timescale shorter than 
Ta.crit ~ 810 yrs is needed in order to pass through the 2:1 
resonance. For example, if we assume that the inner planet 
migrates on a timescale of 2000 years, the outer planet has 
to migrate with a timescale 



^a.2,crit 

< 576 yrs. 



(3) 



We have run several iV-body simulations to explore the 
large parameter space and confirm the above estima te. The 
code used is similar to that presented in Rein fc Papaloizoul 
(2009) and uses a fifth order Runge-Kutta as well as a 
Burlish Stoer integrator, both with adaptive time-stepping. 
Different modules deal with migration and stochastic forc- 
ing. Non conservative forc es are c alculated according to 
the procedure presented in iLee fc P eale (2002) where the 
migration and eccentricity damping timescales = |a/a| 
and Te — |e/e| are imposed for each planet individually. 
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Stochastic forcing was employe d only for one run i n this 
paper and we refer the reader to lRein fc Papaloizoul ()2009[ ) 
for more details on the implementation. 

We place both planets on circular orbits at ai = 1 AU 
and a2 = 2 AU initially. The migration timescale for the 
inner planet is fixed at Ta.i — 2000 yrs, while the migration 
timescale for the outer planet is varied. In figure [1] we plot 
the period ratio P2 /Pi as a function of time t for different 
migration timescales Ta^2- There is a sharp transition of 
the final resonant configuration from 2:1 to 3:2 at around 
Ta,2 ~ 565 yrs. Note that this value agrees extremely well 
with the analytic estimate. 

The above results show that if the planets begin with 
a period ratio exceeding two, in order to get them into the 
observed 3:2 resonance, the relative migration time has to 
be shorter than what is obtained from the standard theory 
of type II m igration applied to these planets in a standard 
model disc ([Nelson et all l2000f ). In that case one expects 
this timescale to be ~ 10^ yrs. However, it is possible to 
obtain the required shorter migration timescales in a mas- 
sive disc in which the pla nets migrate in a type III regim e 
(see e.g. iMasset fc Papaloizoull2003l : IPeplihski et al.ll2008^ . 
In this regime, the surface density distribution in the co- 
orbital region is asymmetric, leading to a large torque which 
is able to cause the planet to fall inwards on a timescale 
much shorter than the disc evolution time obtained for type 
II migration. We explore this possibility in more detail in 
the next section. 



3. Hydrodynamical Simulations 

Two dimensional, grid based hydrodynamic simulations of 
two gravitationally interacting planets which also undergo 
interaction with an accretion disc have been performed in 
order to test the rapid migration hypothesis. The simula- 
tions performed here are similar in concept to those per- 
formed by [Snellgrove et al. (2001) of the resonant coupling 
in the GJ876 system that may have been induced by orbital 
migration resulting from interaction with the protoplane- 
tary disc. We hav e performed studies using the FARGO 
code (|Massetl2000( l with a modified locally isothermal equa- 
tion of state. Those runs are indicated by the letter F. 

We have also run simulations including viscous heating 
an d radiative tra i isport using the RH2D code, very similar 
to iKlev fc Cridal (|2008D . Those runs are indicated by the 
letter R. 



3.1. Equation of state 

The outer planet is likely to undergo rapid type III migra- 
tion and the co-orbital reg ion will be very asymm etric. We 
found, in accordance with IPeplihski et all (|2008l ). that the 
standard softening description does not lead t o convergent 
resu lts in this case (for a comparative study see lCrida et al.l 
l200a ). Due to the massive disc, a large density spike devel- 
ops near the planet. Any small asymmetry will then gen- 
erate a large torq ue, leading to erratic r esults. We follow 
the prescription of IPeplihski et al.l (|2008D and increase the 
sound speed near the outer planet since the locally isother- 
mal model breaks down in the circumplanetary disc. The 
new sound speed is given by 



where Ts and rp are the distance to the star and the outer 
planet, respectively, fls and f2p are the Keplerian angular 
velocity and hg and hp are the aspect ratio, both of the 
circumstellar and circumplanetary disc, respectively. The 
parameter n is chosen to be 3.5 and the aspect ratio of 
the circumplanetary disc is hp ~ 0.4. The sound speed has 
not been changed in the vicinity of the inner, less massive 
planet as it will not undergo type III migration and the 
density peak near the planet is much smaller. 

3.2. Radiation hydrodynamic simulations 

In the radiative simulations we go beyond the locally 
isothermal approximation and include the full thermal en- 
ergy equation, which takes into account the generated vis- 
cous heat and radiative transport 



dt 



-HV- (EcvTu) = -pV-u-|-D-Q-2i?V-F, (5) 



{hsTs 



{hpTpY 



-l/n 



(4) 



where u = {ur,u^) is the two-dimensional velocity, S the 
surface density, p the (vertically averaged) pressure, T the 
(mid-plane) temperature of the disc, and Cv is the specific 
heat at constant volume. On the right hand side, the first 
term describes compressional heating, D the (vertically av- 
eraged) dissipation function, Q the local radiative cooling 
from the two surfaces of the disc (oc T^^), and F denotes the 
two-dimensional radiative flux in the (r, (^)-plane. We cal- 
culate the effective temperature of the disc by T^g — tT"^, 
where r is the vertical optical depth of the disk. To calcu- 
late r, we use a vertically averaged opacity. 

For numerical stability we integrate the diffusion part 
implicitly. Furthe r details abou t solvi ng the energy equa- 
tion are stated in lKlev fc Cridal (l2008t). an d the code RH2D 
is described in more detail in 'iKle^ (Il999h . In this radia- 
tive formulation the sound speed is a direct outcome of the 
simulations and we do not use equation 3] 

3.3. Initial configuration and computational set up 

We use a system of units in which the unit of mass is the 
central mass A/,, the unit of distance is the initial semi- 
major axis of the inner planet, ri, and the unit of time is 
(GAf^/rJ)"-"^/^, thus the orbital period of the initial orbit 
of the inner planet is 27r in these dimensionless units. 

The parameters for some of the simulations we have 
conducted, as well as their outcomes are given in table [H 
the caption being self-explanatory. 

In all simulations a smoothing length b — 0.6 H has 
been adopted, where H is the disc thickness at the planet 
position. This corresponds to about 4 zone widths in a simu- 
lation with resolution of 768 x 768. The role of the softening 
parameter b that is used in two dimensional calculations is 
to take account of the smoothing that would result from the 
vertical structure of the disc in three dimensional calcula- 
tions (e.g. Masset et al 2006). We found that the migration 
rate of the outer planet is independent of the smoothing 
length if the sound speed is given by equation 3) 

The planets are initialised on circular orbits. The mass 
ratios adopted are those estimated for HD45364. We as- 
sume there is no mass accretion so that these remain fixed 
in the simulation (see discussion in iPeplinski et al. 200j, 
and also section [4]). The initial surface density profile is con- 
stant and tests indicate that varying the initial surface den- 
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Table 1. Parameters for some of the hydrodynamic simulations: The first column labels the run, the second and the 
third columns give the mass ratios relative to the central mass of the inner and outer planets, respectively. The fourth 
column gives the initial ratio of the semi-major axes of the two planets. The fifth column gives the initial disc aspect 
ratio, the sixth the softening length in units of the local scale height, the seventh the constant dimensionless viscosity 
and the eighth column gives the initial surface density. The ninth and tenth columns give the inner and outer boundary 
of the 2D grid, respectively. The eleventh and twelfth columns give the radial and azimuthal grid resolutions used and 
finally the simulation outcome is indicated in the final column (see text). 



sity profile does not change the outcome very much. Non- 
reflecting boundary conditions have been used throughout 
this paper. 

3.4. Simulation results 

We found several possible simulation outcomes. These are 
indicated in the final column of table [T] which gives the res- 
onance obtained in each case with a final letter D denoting 
that the migration became ultimately divergent resulting 
in the loss of commensurability. 

We found that convergent migration could lead to a 
2:1 resonance that was set up in the initial stages of the 
simulation (F3). Cases providing more rapid and consistent 
convergent migration than F3 could attain a 3:2 resonance 
directly (F1,F2,F4,F5). 

Thus avoidance of the attainment of sustained 2:1 com- 
mensurability and the effective attainment of 3:2 commen- 
surability required a rapid convergent migration, as pre- 
dicted by the N-body simulations and the analytic esti- 
mate (see section [2?2)) . That is apparently helped initially 
by a rapid inward migration phase of the outer planet that 
approaches a type III migration. The outer planet went 
through that phase in all simulations with a surface density 
larger than 0.00025 and therefore a 3:2 commensurability 
was obtained. Thus, a surface density compa rable to the 
minimum solar nebula (MMSN. lHavashilll98lD . as used in 
simulation F3, is not sufficient to allow the planets undergo 
rapid enough type III migration. As soon as the planets ap- 
proach the 3:2 commensurability, type III migration stops 
due to the interaction with the inner planet and the outer 
planet starts to migrate in a standard type II regime. 

This imposes another constraint on the long term sus- 
tainability of the resonance. The inner planet remains em- 
bedded in the disc and thus potentially undergoes a fairly 
rapid inward type I migration. If the type II migration rate 
of the outer planet is slower than the type I migration rate 
of the inner planet, then the planets diverge and the reso- 
nance is not sustained. Note that a precise estimate of the 
migration rates is not possible at late stages, as the planets 
interact strongly with the density structure imposed on the 
disc by each other. 

Accordingly, outflow boundary conditions at the inner 
boundary that prevent the build up of an inner disc are 
more favourable to the maintenance of a 3:2 commensura- 
bility because the type I migration rate scales linearly with 



the disc surface density. However, those are not presented 
here, as the effect is small and we stop the simulation before 
a large inner disc can build up near the boundary. 

The migration rate for the inner planet depends on the 
aspe ct ratio h. It is dec reased for an increased disc thick- 
ness ()Tanaka et al.|[2002h . That explains why models with a 
large disc thickness (F5) tend to stay in resonance whereas 
models with a smaller disc thickness (F2 , F4) tend towards 
divergent migration at late times. The larger thickness is 
consistent with radiative runs (see below). 

As an illustration of the evolution from a typical config- 
uration (F5) that forms and maintains a 3:2 commensura- 
bility during which the orbital radii contract by a factor of 
at least ~ 2 we plot the evolutions of the semi-major axes, 
the period ratio and eccentricities in figure [2] and provide 
a surface density contour plot in figure [3] after 100 initial 
inner planet orbits at which the 3:2 commensurability has 
just been established. 

The eccentricity peak in figure [2] at t 300 is due to 
passing through the 2:1 commensurability. At t ~ 800 the 
3:2 commensurability is reached and maintained until the 
end of the simulation. The surface density in this simulation 
is app roximately 5 ti mes higher compared to the MMSN at 
1 AU (lHavashilll98l . 

We further present an illustration of the evolution of 
the semi major axes and eccentricities, as well as a surface 
density plot from a simulation (F4) that does form a 3:2 
commensurability, but loses it because the inner planet is 
migrating too fast in figures S] and One can see that a 
massive inner disc has been piled up. This, and and the 
small aspect ratio oi h = 0.04 makes the inner planet go 
faster than the outer planet which has opened a clear gap. 
The commensurability is lost at t ~ 1200. 

The radiative run Rl also shows capture in a 3:2 res- 
onance which is sustained. Before embedding the planets 
the disc is first brought into radiative equilibrium where 
viscous heating is balanced by radiative cooling. During 
this evolution the surface density is kept constant as in 
the local isothermal models. The temperature of the disc 
is obtained directly from the model, and the equilibrium is 
equivalent to approximately H/r = 0.075 at one AU. The 
planets are then embedded at the same initial positions 
as in the other models (F1-F5). We plot the evolution of 
the semi major axes, period ratio and eccentricities in fig- 
ure [6l When reaching the location of the 2:1 resonance the 
eccentricities of both planets increase initially but decline 
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Fig. 2. The semi-major axes (top), period ratio P2/P1 
(middle) and eccentricities (bottom) of the two planets 
plotted as a function of time in dimensionless units for run 
F5 with disc aspect ratio oi h = 0.07. In the bottom panel 
the upper curve corresponds to the inner planet. 
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Fig. 3. A surface density contour plot for simulation F5 
after 100 orbits at the end of the type III migration phase. 
The outer planet establishes a definite gap while the inner 
planet remains embedded at the edge of the outer planet's 
gap. 
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Fig. 4. The semi-major axes (top), period ratio P2/P1 
(middle) and eccentricities (bottom) of the two planets 
plotted as a function of time in dimensionless units for run 
F4 with a disc aspect ratio oih = 0.04. In the bottom panel 
the upper curve corresponds to the inner planet. 
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Fig. 5. A surface density contour plot for simulation F4 af- 
ter 150 orbits after the planets went into divergent migra- 
tion. The inner planet is embedded and interacts strongly 
with the inner disc. The simulation uses a ID grid for 
0.04 <r < 0.25. 



subsequently after the resonance is transversed. Capture 
into 3:2 occurs smoothly during the inward migration of 
the planets. The final eccentricities reached in this model 
are around ei ~ 0.03 and 62 — 0.01. 

The initial type III migration rate is slower when com- 
pared to run F5 because the surface density is lower. Note 
that after the the 2:1 resonance is passed, the outer planet 
migrates in a type II regime and thus, the capture into 
the 3:2 resonance appears later. However, the orbital pa- 
rameters measured at the end of the simulations are very 
similar to any other run that we have performed. This indi- 
cates that the parameter space that is populated by these 
kind of planet disc simulation is very generic. 



We tested the effect of disc dispersal at the late stages 
in our models to evolve the system self-consistently to 
the present day. In model F5 after t = 2000 we allow 
the disc mass to exponentially decay on a timescale of 
Tdis ^ 2000. Note that this ti mescale is shorter tha n 
the photo-evaporation timescale ([Alexander et al.l |2006[ ). 
However, this scenario is expected to give a stronger effect 
than a long timescale (jSandor fc Kiev [20061 ). In agreement 
with those authors, we found that the dynamical state of 
the system does not change for the above parameters. At a 
late stage, the resonance is well established and the planets 
undergo a slow inward migration. Strong effects are only 
expected if the disc dispersal happens during the short pe- 
riod of rapid type III migration, which is very unlikely. We 
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Fig. 6. The semi- major axes (top), period ratio P2/P1 
(middle) and eccentricities (bottom) of the two planets 
plotted as a function of time in dimensionless units for the 
radiative run Rl. 

observed that the eccentricities show a trend to decrease 
and the libration amplitudes to slightly increase during the 
dispersal phase. However, these changes are not different 
than what has been observed in runs without a disappear- 
ing disc. 

4. Other scenarios for the origin of HD45364 

In the above discussion we have considered the situation 
when the planets attain their final masses while having a 
wider separation than that required for a 2:1 commensura- 
bility and found that convergent migration scenarios can be 
found that bring them into the observed 3:2 commensura- 
bility by disc planet interaction. However, it is possible that 
they could be brought to their current configuration in a 
number of different ways as considered below. Nonetheless 
it is important to note that because the final commensu- 
rable state results from disc planet interactions it should 
have similar properties as those described above when mak- 
ing comparisons with observation. 

It is possible that the solid cores of, either both planets, 
or just the outer planet approached the inner planet more 
closely than the 2:1 commensurability before entering the 
rapid gas accretion phase and attaining their final masses 
prior to entering the 3:2 commensurability. Although it is 
difficult to rule out such possibilities entirely, we note that 
the cores would be expected to be in the super earth mass 
range where in general closer commensurabilities than 2:1 
and even 3:2 are found for typical type I migration rates 
(e.g. iPapaloizou fc Szuszkiewicd [20051 : ICresswell fc Nelsonl 
I2008D . One may also envisage the possibility that the solid 
cores grew in situ in a 3:2 commensurability, but this would 
have to survive expected strongly varying migration rates 
as a result of disc planet interactions as the planets grew 
in mass. 

An issue is whether the embedded inner planet is in 
a rapid accretion phase. The onset of the rapid accretion 
phase (also called phase 3) o ccurs when the core and enve- 
lope mass are about equal (jPollack et al.|[l996l ). The to- 



tal planet mass depends at this stage on the boundary 
conditions, here determined by the circumplanetary fiow. 
When these allow the planet to have a significant convec- 
tive envelope the transition to rapid accre tion may not oc - 
cur until the planet mass exceeds 6OM0 (jWuchterllll993[ ). 
which is the mass o f the inner plan et (see also model J3 of 
Pollack et al.lll99^ and models of IPapaloizou fc TerquemI 
19991) . Because of the above results, it is not unreasonable 
that the inner planet is not in a rapid accretion phase. 

Finally we remark that protoplanetary discs are be- 
lieved to maintain turbulence i n some part of thei r struc - 
ture. Using the prescription of iRein k. Papaloizoul (|2009[ ). 
we can simulate the turbulent behaviour of the disc by 
adding additional forces to an N-body simulation. These 
forces will ultimately eject the planets from the 2:1 reso- 
nance should that form. Provided they are strong enough 
this can happen within the lifetime of the disc, thus mak- 
ing a subsequent capture into the observed 3:2 resonance 
possible. 

We have confirmed numerically that such cases can oc- 
cu r for moderately large dif fusion coefficients (as estimated 
bv lRein fc Papaloizoull2009 l . However, this outcome seems 
to be the exception rather than the rule. Should the 2:1 
resonance be broken, a planet-planet scattering event ap- 
pears to be more likely. In all the simulations we performed, 
we found that only a small fraction of systems (1% - 5%) 
eventually end up in a 3:2 resonance. 

5. Comparison with observations 

Table [2] lists the orbital parameters that ICorreia et al.l 
(|2009D obtained from their best statistical fit of two 
Keplerian orbits to the radial velocity data. The radial ve- 
locity data collected so far is insufficient to detect the pres- 
ence of any interactions between the planets. This is borne 
out by the fact that the authors did not obtain an improve- 
ment in terms of minimising 'x^ when a 3-body Newtonian 
fit rather than a Keplerian fit was carried out. However, a 
stability analysis supports the viability of the determined 
parameters as these lie inside a stable region of orbital pa- 
rameter space. The best fit solution shows a 3:2 mean mo- 
tion commensurability although no planet-planet interac- 
tions have been observed from which this could be inferred 
directly. It is important to note that the large minimum 

associated with the best fit that could be obtained in- 
dicates large uncertainties which are not accounted for by 
the magnitudes of the errors quoted in their paper. 

There are two main reasons why we believe that the 
effective errors associated with the observations are larger: 
First, many data points are clustered and clearly do not 
provide a random time sampling. This could result in cor- 
relations which effectively reduce the total number of in- 
dependent measurements. Second, the large minimum 
value being three times the quoted observational error in- 
dicates that there are additional effects (e.g. sunspots, addi- 
tional planets, etc) that produce a jitter in the central star 
(M. Mayor, private communication) which then enhances 
the effective observational error. 

No matter what process generates the additional noise, 
we can assume in a first approximation that it follows a 
Gaussian distribution. We then have to conclude that the 
effective error of each observation is a factor ^ 3 larger than 
reported in order to account for the minimum given the 
quoted number of ~ 60 independent observations. Under 
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Uorreia et al. ( zUUa ) 


Simulation F5 


Parameter 


Unit 


D C 


D C 


Af sin i 




U.io</ U.DOiy 


n 1 Qvo o dK^Ci 
U.io/2 U.DOiy 


M, 




0.82 


0.82 


a 


[AU] 


0.6813 0.8972 


0.6804 0.8994 


e 


0.17 ±0.02 0.097 ±0.012 


0.036 0.017 


A 


[deg] 


105.8 ± 1.4 269.5 ±0.6 


352.5 153.9 




[deg] 


162.6 ±6.3 7.4 ±4.3 


87.9 292.2 






2.79 


2.76* (3.51) 


Date 


[JD] 


2453500 


2453500 



Table 2. Orbital p arameters of HD45364b and HD45364c. The observed values including error-bars are taken from 
ICorreia et al.l (|2009t ). Note that vu is not well constrained for nearly circular orbits. The radial velocity amplitude of the 
simulation marked with * has been scaled down by 8%. 




Observations 



Correia et al 




Simulation F5 
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53200 



53400 



53600 53800 
JD-2400000 [days] 
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54200 
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54600 



Fig. 7. Comparison of different orbital solutions. See text for a description of the diffe rent models. Radial velocity 
measurements (red points) and the published orbital solution (green curve) are taken from ICorreia et al.l (I2009D . 



these circumstances we would then conclude that any fit 
with \/x^ in the range 2 — 4 would be an equally valid 
possibility. 

In the following, we show that indeed a variety of or- 
bital solutions match the observed RV data with no statis- 
tically sigi iificant difference wh en compared to the quoted 
best fit bv ICorreia et al.l (|2009t ). As one illustration, we use 
the simulation F5 obtained in section |3l We take the or- 
bital parameters at a time when the orbital period of the 
outer planet is closest to the observed value and integrate 
them for several orbits with our N-body code. The solu- 
tion is stable for at least one million years. There are only 
two free parameters available to fit the refiex motion of the 
central star to the observed radial velocity: the origin of 
time (epoch) and the angle between the line of sight and 
the pericentre of the planets. 

We can assume that the planet masses and periods are 
measured with high accuracy and that only the shape of 
the orbit contains large errors. Our best fit results in an 
unreduced \f)^ value of 3.51. According to the above dis- 
cussion, this solution is statistically indistinguishable from 
a solution with ~ 2.8. 



It is possible to reduce the value of even further, 
when assuming that the planet masses are not fixed. In 
that case, we have found a fit corresponding to \/)^ — 2.76 
where we have reduced the radial velocity amplitude by 8% 
(effectively adding one free parameter to the fit). This could 
be explained either by a heavier star, less massive planets 
or a less inclined orbit. The orientation of the orbital plane 
has been kept fixed {i = 90°) in all fits. 

The results are shown in figure [71 The blue curve cor- 
responds to the outcome of simulation F5 and we list the 
or bital parameters in ta blejH We also plot the best fit given 
bv ICorreia et al.l (|2009i) for comparison (green curve). In 
accordance with the discussion given above it is very diffi- 
cult to see any differences in the quality of these fits which 
indeed suggests that the models are statistically indistin- 
guishable. 

However, there is an important differenc e between the 
fits ob tained from our simulations and that of lCorreia et al.l 
(|2009D . Our models consistently predict lower values for 
both eccentricities ei and e2 (see tabled]). Furthermore the 
ratio of eccentricities exjei is higher than the previously 
reported value of 1.73. Note however, that the eccentricities 
are oscillating and the ratio is on average ~ 3. 
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reported fit by ICorreia et al.l (120091 ) and the new fit ob- 
tained from simulation F5. 
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Fig. 9. Evolution of the resonant angles 0i = 2Ai — 3A2 + 
^^ii 02 = 2Ai — 3A2 + •072, ci7i — 1372 for simulation F5 using 
the parameters given in table [21 The angles are measured 
in radians. 

This in turn results in a different li bration pattern: the 
slow libration mode (see figure 3a in ICorreia et all |2009[ ) 
that is associated with oscillations of the angle between the 
two apsidal lines is absent, as shown in figure [9l Thus there 
is a marked difference in the form the interaction between 
the two planets takes in the case of this fit. 

It is hoped that future observations will be able to re- 
solve this issue. The evolution of the difference in radial 
velocity over the next few yea rs between our fit fro m simu- 
lation F5 and the fit found by ICorreia et al.l (|2009 f) is plot- 
ted in figure [Hj There is approximately one window each 
year that will allow further observations to distinguish be- 
tween the two models. The large difference at these dates is 
not due to the secular evolution of the system, but simply 
a consequence of smaller eccentricities. 

All runs that yield a final configuration with a 3:2 MMR 
are very similar in their dynamical behaviour, as found 
for model F5. They are characterised by an antisymmet- 



ric (1372 — wi ^ n) state with relatively small libration of (jji 
and 02 • This is a generic outcome of convergent disk-planet 
migration scenarios. 

6. Conclusions 

The planets in the multi-planetary system HD45364 are 
most likely in a 3:2 mean motion resonance. This poses 
interesting questions on it's formation history. Assuming 
that the planets form far apart from each other and migrate 
with moderate migration rate, as predicted by standard 
planet formation and migration theories, the most likely 
outcome is a 2:1 mean motion resonance, contrary to the 
observation of a 3:2 MMR. 

In this work, we investigated a possible way around this 
problem by letting the outer planet undergo a rapid inward 
type III migration. We presented an analytical estimate and 
performed N-body as well as hydrodynamical simulations. 

Wc found that it is indeed possible to form a 3:2 
MMR and avoid the 2:1 resonance, thus resembling the 
observed planetary system using reasonable disc parame- 
ters. Hydrodynamical simulations suggest that the system 
is more likely to sustain the resonance for large aspect ra- 
tios, as the migration of the inner planet is slowed down, 
thus avoiding divergent migration. 

Finally, we used the orbital configuration found in the 
hydrodynamical formation scenario to calculate a radial ve- 
locity curve. This curve has then been compared to obser- 
vations and the resulting fit has an identical value than 
the previously reported best fit. 

Our solution is stable for at least a million years. It is 
in a dynamically different state, both planets having lower 
eccentricities and a different libration pattern. This is the 
first time that planet migration theory can predict a precise 
orbital configuration of a multiplanetary system. 

This might also be the first direct evidence for type III 
migration if this scenario turns out to be true. The system 
HD45364 remains an interesting object for observers, as the 
differences between the two solutions can be measured in 
radial velocity within a couple of years. 
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